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Abstract 

Equations for the electron Green’s function of the two-dimensional 
Hubbard model, derived using the strong coupling diagram technique, are 
self-consistently solved for different electron concentrations n and tight- 
binding dispersions. Comparison of spectral functions calculated for the 
ratio of Hubbard repulsion to the nearest neighbor hopping U/t = 8 with 
Monte Carlo data shows not only qualitative, but in some cases quanti¬ 
tative agreement in position of maxima. General spectral shapes, their 
evolution with momentum and filling in the wide range 0.7 < n < 1 are 
also similar. At half-filling and for the next nearest neighbor hopping 
constant t' = —0.3t the Mott transition occurs at Uc ~ 7A/8, where A is 
the initial bandwidth. This value is close to those obtained in the cases of 
the semi-elliptical density of states and for t' — 0. In the case U = 8t and 
t' = —0.3t the Mott gap reaches maximum width at n = 1.04, and it is 
larger than that at t' = 0 for half-filling. In all considered cases positions 
of spectral maxima are close to those in the Hubbard-I approximation. 

1 Introduction 

The fermionic Hubbard model is one of the main models used for the description 
of strong electron correlations in crystals such as cuprate perovskites and heavy 
fermion compounds. The strong coupling diagram technique [HI El HIS] is one 
of the approximate approaches used for investigating the model. The method 
is based on the serial expansion in powers of the electron kinetic energy. The 
elements of the arising diagram technique are on-site cumulants of electron 
creation and annihilation operators and hopping lines connecting cumulants on 
different sites. As in the diagram technique with the expansion in powers of an 
interaction, in the considered approach the linked-cluster theorem allows one to 
discard disconnected diagrams and to carry out partial summations in remaining 
connected diagrams. As a consequence the electron Green’s function is expressed 
in the form of the Larkin equation containing the initial electron dispersion and 
the irreducible part - the sum of all irreducible diagrams without external ends. 
In spite of the fact that by the construction the approach is intended for the 
case of strong coupling, when the Hubbard repulsion U is approximately equal 
to or larger than the initial bandwidth A, it gives the correct result in the limit 


1 


[/ —>■ 0. Hence the approach provides an interpolation between the limits of 
weak and strong correlations. 

Using this method, in recent work [5] equations for the electron Green’s 
function were obtained by keeping terms of the lowest two orders in the irre¬ 
ducible part of the Larkin equation. Self-consistent calculations performed for 
the semi-elliptical density of states (DOS) showed that at half-filling the ap¬ 
proximation describes the Mott transition, which occurs at Uc = •\/3A/2. This 
value coincides with the critical repulsion obtained for the same DOS in the 
Hubbard-III approximation [7] . In [S] the same method was applied to the half- 
filled two-dimensional (2D) Hubbard model with nearest neighbor form of the 
kinetic energy (the t-U model). In this case the Mott transition is observed 
at Uc ~ 7A/8 - the value, which is very close to that for the semi-elliptical 
DOS. The used approach allows one to describe not only the Mott transition, 
but also another important feature of the model ~ spectral continua, which lead 
to finite widths of spectral maxima. Spectral functions, DOS and momentum 
distributions calculated in the t-U model were in qualitative and in some cases 
quantitative agreement with Monte Carlo data. 

In the present work the equations derived in are self-consistently solved for 
different electron concentration n and for the kinetic energy containing nearest 
and next nearest neighbor hopping terms (the t-t'-U model). Similar dispersions 
are frequently used for the description of hole-doped cuprates with the ratio of 
the next nearest neighbor to nearest neighbor hopping constants close to t'/t = 
—0.3 (see, e.g., i)- At half-filling and for this ratio of t' jt the Mott transition is 
observed again at Uc ~ 7A/8, close to the values mentioned above. Comparison 
of calculated spectral functions with available Monte Carlo results for the case 
U = 8t and t' = 0 shows not only qualitative, but in some cases quantitative 
agreement in positions of maxima. General spectral shapes, their evolution 
with wave vector and concentration in the wide range 0.7 < n < 1 are also 
close. These facts and relative simplicity of calculations give promise that the 
considered approach can be applied for investigating multi-band generalizations 
of the Hubbard model, for example, such as used for the description of iron 
pnictides or transition-metal oxides. It was found that the behaviour of the Mott 
gap with changing n is different in the t-U and t-t'-U models. In all considered 
cases positions of maxima in spectral functions are close to the frequencies of 
d-function peaks in the Hubbard-I approximation m- 

2 Main formulas 

In this section some equations derived in are reproduced and converted to 
the form, which is convenient for calculations. The electron Green’s function 


G(n'T',lr) = (Ta„/„(r')a„^(T)) 


( 1 ) 


is considered, where the angular brackets denote the statistical averaging with 
the Hubbard Hamiltonian 



tnn' is the hopping constants, the operator aj,^ creates an electron on the site n 
of the 2D square lattice with the spin projection a = ±1, the electron number 
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operator rina = Oncr^ncr, M is the chemical potential, T is the time-ordering op¬ 
erator which arranges operators from right to left in ascending order of times r, 
o.n( 7 {T) = exp(i?T)anCT exp(—ifr) anddnCT(T) = exp(iJr)aJjg. exp(—iJr). Green’s 
function o does not depend on the spin projection, and it was omitted in the 
function notation. 

In the strong coupling diagram technique, after the summation of all dia¬ 
grams the Fourier transform of Green’s function ([T]) acquires the form of the 
Larkin equation 


G(k,^a;^) 


K{'k,iuji) 

1 -ti^K{k,iuJi)' 


(3) 


where k is the 2D wave vector, w; = {21 + 1)ttT is the Matsubara frequency with 
the temperature T, tk = ^jjexp[ik(n — n')]tnn/ and K{k,iuji,) is the irreducible 
part - the sum of all irreducible diagrams without external ends. A diagram is 
said to be an irreducible one if it cannot be divided into two disconnected parts 
by cutting some hopping line inn'- 

In this work K (k, iuji) is approximated by the sum of two terms of the lowest 
orders in powers of tk- These terms contain cumulants of the first and second 
orders. The term with the second-order cumulant includes also a hopping-line 
loop. Using the possibility of the partial summation of diagrams we transform 
this bare hopping line into the dressed one. 


0{nT, n'r') = tnn'S{T - t') + ^ tnm 

mm' 



where /3 = 1/T. With this substitution, for a chemical potential in the range 


/r » T, [7 — ^ » T (5) 

the expression for the K{'k,iuji,) reads 

K{iuJi) = i [goiiiuji) + gi2{iuJi)] 

+ ^F'^{iuJi)(j){iuJi) - ^F{iuji) - ^J{iuJi), ( 6 ) 

where 


5 oi(jw;) = {iwi -I- , gi2{iuii) = {itoi + fJ-- Uy^, 

F{ioJi) = goiiiuJi) - 512 ( 10 ;;), 

J{iuji) = goiiioJi) - F{iuJi)gi2{iuJi), (7) 

si=T '^J{iuji)(j){iuji), s 2 =t'^ F{iu}i)4>{iu}i), 

i i 

(j){iuji) = N~'^ ^ tlG{k, iuJi) 

k 

and N is the number of sites. We chose tnn — 0 and used this relation in deriv¬ 
ing dni). Equations ©, © and © form the closed set of equation for calculating 
Green’s function Notice that irreducible part © does not depend on mo¬ 
mentum as the approximation does not take into account the interactions of 
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electrons with the magnetic ordering, spin and charge fluctuations. In the used 
approach these interactions are described by the sum of diagrams containing 
ladders, which correspond to the dynamic spin and charge susceptibilities m- 
For further consideration we perform the analytic continuation to real fre¬ 
quencies w and transform the sums si and S 2 to the form more convenient for 
calculations using the spectral representations, Poisson summation formulas and 
inequalities ®, 


Si 


S2 


2U — oj — ^ , .dui 

u {uj + y, — u y TT 

, . doj 
Im0(a;) —, 

TT 


Jo uiuj + fj.y 

Im(^(a;) du 

-c^UJ + ^i-U TT 


duj 

UJ + jj, TT 


From and © we get 


( 8 ) 


4>{uj) = {oj yf{uT + ^-U)‘^N ^ 


X s (w -i- y)^{u) ^ 


—tk + ^ + ~ 1) ) (^ + /^) 

X (w -|- /I — [/) —— I (w -|- /I — C/)^ 


+U (w -I- ^) ) -I- 


-1 


(9) 


A self-consistent solution (j){uj) of ® and ® is used in ® and ® for calculating 
Green’s function. 

As follows from frequencies ui = and U — ^ are somehow singled 

out. At half-filling and in the case of a DOS, which is symmetric with respect 
to fi = U/2, the sums si and S 2 vanish. As follows from ®, in this case 4>{uj) 
behaves near uj = ±U/2 as (see also [S]) 



Substituting this result into ® we see that Im Ar(a;) diverges as (w ± 17/2)“^. 
This fact presents a problem because up to a multiplier ImAr(a;) is a spectral 
function for momenta meeting the condition tk = 0 [see ®], which due to the 
divergence does not satisfy neither the normalization condition 


Im K{u}) duj = —TT 


nor the Kramers-Kronig relation 
ReA:(a;) 


ImAr(a;') duj 


OJ' — OJ 


( 10 ) 


( 11 ) 
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Figure 1: Spectral functions in the t-U model for U = 8t,n = 0.93, a 6 x 6 lattice 
and momenta shown in the panels. Blue solid lines are our results, red dashed 
lines are Monte Carlo simulations for T = 0.5t [12], arrows indicate locations of 
(^-function peaks of the Hubbard-I approximation. 

To overcome this difficulty in [5] we introduced an artificial broadening. How¬ 
ever, if an initial DOS is not symmetric (as for the case t' ^ 0) and/or n 7 ^ 1 
the divergence disappears. Indeed, in this case si and S 2 in (jHj) are nonzero, and 
from Q we get near w = —/i and U — 

2 

(j){oj fv -fi) fv —{ui + nf, 

S2 

cj}(u} KiU — fi) Ki —(o; -f /J. — UY. 

S2 

Here (j){oj) is purely real and, therefore, IrccK = 0. This poses another difficulty: 
as follows from calculations, for S 2 7 ^ 0 narrow gaps with nearly vertical walls 
appear in ImiF(w) around ui = —fi and U — fi. Thanks to their narrowness and 
shape, these unphysical gaps can be easily eliminated by a linear interpolation 
between tops of the walls. 

3 Results and discussion 

Results of this section were obtained by the self-consistent solution of (Hj) and (l9|) 
for a given /r in the t-U or t-t'-U models. The solution was found by iteration. 
Such obtained /'(w), si and S 2 were used for calculating ImiF(w) from ([S]). 
After removing the gaps near lo = —/i and {7 — /r with the interpolation and 
normalization (nni, lmK{uj) is used for calculating YieK{uj) from (fTTll . which 
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Figure 2: Same as in Fig. [U but for n = 0.7 and a 8x8 lattice. The Monte 
Carlo data for T = 0.5t are taken from |13| . 


allows one to calculate Green’s function m for arbitrary k. Such obtained 
spectral functions A{kui) = —7r“^Im G'(k(u;) and DOS p{ui) = A(ka;) 

satisfy the usual normalization conditions with good accuracy. 

To check the validity of the used approach its results were compared with 
data of Monte Carlo simulations carried out in [12] and [13] . This comparison 
is shown in Figs. [T] and [2] Here and below t and the intersite distance are 
set as units of energy and length, respectively. The electron concentration was 
determined from the relation 


n = 


H(ka;) 

N ^ j_^ exp(a;/T) + 1 


where T was set equal to 0, since in the used approximation equations do not 
depend on T and actually correspond to T = 0. Comparing our and Monte Carlo 
results it should be borne in mind that due to the sign problem the latter were 
obtained at comparatively high temperature T = 0.5t. This leads to increased 
widths of maxima in comparison with zero-temperature results. We did not 
try to fit widths of maxima in our spectra to Monte Carlo data by an artificial 
broadening. This was done to demonstrate real spectral shapes, which appear 
in the present approach. Besides, Monte Carlo spectra for real frequencies 
were obtained by analytic continuation of imaginary-frequency data using the 
maximum-entropy method. This method introduces additional inaccuracies in 
spectral shapes. 

Keeping these remarks in mind, from Figs. |T] and [5] we can conclude that 
the used approach gives spectral functions in qualitatively and in some cases 
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Figure 3: Densities of states in the t-U model for U = 8t, n = 1 (blue solid 
line), 0.99 (red dashed line), 0.94 (olive dash-dotted line) and 0.81 (purple short- 
dashed line). A 240x240 lattice. 

quantitatively agreement with the Monte Carlo data in the wide range of elec¬ 
tron concentrations 0.7 < n < 1 (see also the analogous comparison for the case 
of half-hlling in 0 ). Locations of maxima, general shapes of spectra and their 
variation with k and n are as a rule close. Substantial differences in widths 
of maxima, observed in some spectra, may be related to the above-mentioned 
reasons. A similar picture is observed also in the comparison of Monte Carlo 
m and our results for n = 0.87 (not shown here). 

Arrows in Figs. [1] and [2] indicate frequencies of ^-function peaks in the 
Hubbard-I approximation, 



As for the case of half-filling [5], one can see that in the used approximation 
locations of maxima are close to these frequencies. 

Widths of spectral maxima are essentially increased away from half-filling. 
This may be observed in the evolution of the DOS with changing n in Fig. [21 Due 
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Figure 4: Densities of states in the t-t'-U model for U = 8t, t' = —0.3t, n = 1.26 
(purple dash-dot-dotted line), 1.04 (pink short-dashed line), 1 (blue solid line), 
0.99 (red dashed line) and 0.87 (olive dash-dotted line). A 240x240 lattice. 
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Figure 5: Spectral functions in the t-t'-U model for U = 8t, t' = —0.3t, n = I 
and momenta along the symmetry lines (0,0)-(7r, 7r)-(7r, 0)-(0, 0). A 240x240 
lattice. Red pins on the base lines of each curve indicate locations of 5-function 
peaks in the Hubbard-I approximation. 



0) 

Figure 6: Same as in Fig.jSJ but for n = 1.26. 


to the particle-hole symmetry of the t-U model only cases n < 1 are shown. The 
DOS for n > 1 can be obtained from the respective case n < 1 by the specular 
reflection in the line w = 0. As already noted earlier mm, with deviation from 
half-filling the DOS is redistributed in favor of the Hubbard subband, in which 
the Fermi level is located. The Mott gap decreases monotonously with growing 
|1 — n| and disappears at |1 — n| ~ 0.12. 

Similar behaviour is observed in the DOS of the t-t'-U model shown in Fig. 01 
However, in this case the width of the Mott gap reaches its maximum value at 
n ~ 1.04 rather than at n = 1, as in the t-U model. In the t-t'-U model the 
maximal gap width is somewhat larger - 0.64t at {7 = 8t and t' = —0.3t in 
comparison with 0.54t at t' = 0. The dependence of the gap width on n — 1.04 
is strongly asymmetric in the t-t'-U model. At half-filling, keeping t' = —0.3t 
and decreasing U we found that the gap disappears at Uc~ It = 7A/8. This 
value of the critical repulsion is very close to those found for the semi-elliptical 
DOS [TIE] and in the t-U model [5]. For the used parameters the positions of 
maxima in spectral functions of the t-t'-U model are also close to those in the 
Hubbard-I approximation, as can be seen in Figs. jSjandEl 
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4 Concluding remarks 

In this work, equations for the electron Green’s function of the repulsive Hub¬ 
bard model, derived using the strong coupling diagram technique, were self- 
consistently solved for different electron concentrations n in the two-dimensional 
t-U and t-t'-U models. Terms of the first two orders in the expansion in powers 
oit/U were taken into account in the irreducible part of the equation for Green’s 
function, and the bare internal hopping line in one of these terms was substi¬ 
tuted with the dressed one. Gomparison of spectral functions, calculated in the 
t-U model for the Hubbard repulsion U = 8t, with Monte Garlo data shows not 
only qualitative but in some cases quantitative agreement in positions of max¬ 
ima. General spectral shapes, their evolution with the wave vector and electron 
concentration in the wide range 0.7 < n < 1 are also similar. It was found that 
in the half-filled t-t'-U model with t'/t = —0.3 the Mott transition occurs at 
Uc ~ 7A/8, where A is the initial bandwidth. This critical value is close to 
those found for the semi-elliptical density of states and in the case t' = 0. The 
behaviour of the Mott gap is different in the cases t' = 0 and t' = —0.3t: the 
maximal gap width is larger in the latter case and is reached at n = 1.04 rather 
than at half-filling, as in the case t' = 0. In both cases positions of spectral 
maxima are close to those in the Hubbard-I approximation. 

The conducted comparison with Monte Garlo results and comparative sim¬ 
plicity of calculations give promise that the strong coupling diagram technique 
may be a useful tool in consideration of different generalizations of the Hubbard 
model such as multi-band models used for the description of transition-metal 
oxides or models containing external fields. 
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